Zeros and orthogonality of the Askey- Wilson 
polynomials for q a root of unity^ 

Vyacheslav Spiridonov and Alexei Zhedanov 

Abstract 

We study some properties of the Askey- Wilson polynomials (AWP) when g is a primitive 
. Ath root of unity. For general four-parameter AWP, zeros of the Ath polynomial and the 

I orthogonality measure are found explicitly. Special subclasses of the AWP, e.g., the continuous 

g-Jacobi and big g-Jacobi polynomials, are considered in detail. A set of discrete weight functions 
positive on a real interval is described. Some new trigonometric identities related to the AWP 
are obtained. Normalization conditions of some polynomials are expressed in terms of the Gauss 
sums. 
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§: 1. Introduction 

bX)' The theory of orthogonal polynomials of one variable is based upon the three-term recurrence relation 

13: [8] 

Pn+l{x) +hnPn{x) + UnPn-l{x) = xPn{x), 71 = 1,2,... (1.1) 

■ with the initial conditions 

X'- Po{x) = l, P,ix)=x-bo. (1.2) 

d ' There is a representation theorem concerning orthogonality of the polynomials P„(x) with respect 
to some complex measure fj,{x) in case Un and 6„ are arbitrary complex coefficients; see, e.g., [3, 8], 

Pn{x)Pm{x)d^{x) = hnSnm, (1-3) 

C 

where /i„ = U1U2 ■ ■ - Un and C is some contour in the complex plane. When bn are real and m„ > 0, 
the measure can be chosen real-valued (then C is the real line) and nondecreasing. For such cases it 
can be shown that all zeros Xs, s = 0,1, . . . , N — 1, of any polynomial Pn{x) are real and simple [8]. 
Moreover, these zeros can be used for construction of a discrete orthogonality relation for polynomials 
of the degree lower than N: 

N-l 

Pn{Xs)Pm{Xs)Ws = K^nm, (1-4) 

where 72,, 7n = 0, 1, . . . , — 1 and the weight function is 

PN-l{Xs)P'N{^sy ^ ' 



^Date: April 1996; published in Duke Math. J. 89 (1997), 283-305. 
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The property (1.4) was proven in [4] under the assumption that 6„ are real and Un are positive. 
However, the proof remains vahd even in general cases of complex bn,Un. The only requirement to 
be added is that all zeros Xs must be simple which is not guaranteed if the recurrence coefficients arc 
complex. For example, for fo„ = one finds from the recurrence relation (1.1), P3{x) = x{x'^ — Ui — U2)- 
If U2 — —Ui then all zeros of the polynomial i-*3(x) are located at x — 0. It is thus always necessary 
to check simplicity of the zeros x, in order for (1.4) to be applicable. 

In this paper we present a study of the finite orthogonality property (1.4) of the Askey- Wilson 
polynomials (AWP) for the case when g is a root of unity. Recall that the recurrence coefficients for 
AWP have the form [3, 9] 

Un^^nVn-U K ^ -^n - Vn + a + a~'^ , (1.6) 

where 

_ a(l - g")(l - 6cg"-i)(l - Mg"-i)(l - cdq""-^) 

^" - (l_^g2n-2)(l_^^2n-l) ' ( " ) 

_ a-\l - a6g")(l - acg")(l - arfg")(l - ^g""^) 

- (l-^g2«)(l-^g2n-i) ' (1-8) 

and g = abed. The parameters a, b, c, d are arbitrary complex numbers. The explicit expression of 
AWP in terms of the basic hypergeometric functions is [3, 9] 

where 

n _ a~''{ab;q)n{ac;q)n{ad]q)n . . 

and X = t + 1^^. Here (a; g)„ is the g-shifted factorial defined as (a; g)o = 1 and (a; g)„ = (1 — a)(l — 
aq) ■ ■ ■ (1 — ag"^^) for integer n > 0. The general g-hyper geometric series is defined as follows [9]: 

\ bi,...,bs J ^0 {Q,bu---,bs;q)n 

where the compact notation (ai, a2, • • • , a^; q)n — (oi; q)n • • • {dr] (l)n is used. 

For real q and appropriate restrictions upon the parameters a, 6, c, d the AWP are orthogonal on 
a finite interval of the real axis with a positive continuous weight function [3]. In this case, AWP 
satisfy a difference equation defined upon the hyperbolic grid [3, 15]. A possibility of constructing 
orthogonal polynomials for trigonometric grids has been mentioned briefly in [15]. In this paper we 
concentrate on the case when g is a primitive A'"-th root of unity: 

q = exp{27iiM/N), (1.12) 

where M and are mutually co-prime positive integers, M < N. We consider only finite sets of or- 
thogonal polynomials with discrete measures which are different from the so-called sieved polynomials 
(see, e.g., [1]). 



2. Zeros of the A^th Askey- Wilson polynomial 

In this section we analyze location of zeros for generic A^th order AWP depending on four parameters 
for q a primitive A^th root of unity. 
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Consider the explicit form of P/v(.t) (1.9) for q fixed by (1.12). Obviously we cannot use this 
expression directly, because the factor {q^^ \ (i)n/ (?; q) n in the A^th term of the series is indeterminate. 
Nevertheless, we can apply a limiting procedure taking q = exp(e + 2mM/N) and going to the limit 
e — > 0. Then it is easy to see that 

hm 1^:!^ = lim (-lfq-^(^+m ^ _i. (2.1) 
Using (2.1) one can find that 

Pn(x) ^Dn(i- q)N{at-^;q)N \ _ .2.2) 

^ ^ \ [ah] q)N{ac] q)Niad] q)N ) ' 

Since = in (1.1). we have a finite set of polynomials and it is not necessary to consider 
ambiguities in the higher order terms in the series (1.9). 

The Q'-shifted factorial in (2.2) can be simplified because for any complex parameter a one has 
the identity (which is valid only for q a root of unity) 

(a; g)Ar = 1 - a^. (2.3) 

Indeed, (a; q)^ = (1 — a)(l — aq) ■ ■ ■ (1 — aq^~^) is a A^-th order polynomial in argument a. The zeros 
of this polynomial are 1, g, g^, . . . , whence (2.3) follows. 
We can rewrite (2.2) in the form 

Pjv(x) = i^ + r^-£;jv, (2.4) 

where is a rational function of parameters, 

+ 6^ + + - (a6c)^ - (6crf)^ - (a6d)^ - (acd)^ 

= T^T^k^F ■ 

We call Ej^ a combinatorial invariant because its symmetry with respect to permutation of the 
parameters a, 6, c, d plays an important role in the following. 
From (2.4) it is not difficult to find zeros of Pn{x): 

Xs^rq' + r-^q-', s ^ 0,1, . . . , N - 1, (2.6) 

where for simplicity we assume that r is the root with minimal argument of the equation 

r^ = £;jv/2 + ^£;y4-l, or r"^ = Er,/2 - ^Je^/A - 1. 

Let us point out a curious identity associated with the combinatorial invariant E^: 

(1 - aV)(l - 6^r^)(l - c^r^)(l - d^r^) 



(a^ - r^){b^ - r^)(c^ - r^){d^ - r^) 



= 1. (2.7) 



Now we are ready to formulate the orthogonality relation for AWP. First we need the condition 
of simplicity of the zeros Xg, which is 

En 7^ ±2. (2.8) 
Second, in order for the polynomials Pn{x),0 < n < N, to be well defined we impose the constraints 

g,ah,ac,ad,hc,hd,cd ^ q^ , k — 0,1, . . . , N — 1, (2.9) 
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where q is fixed by (1.12). These restrictions guarantee that Un does not vanish between two zeros 
at n = and n — N or that Un,bn do not blow up (cf. [16]). Under the conditions (2.9) one can 
calculate 

N _ -N 

P'A^s) = N— (2.10) 

Using the same e — > limiting procedure as in (2.1) for (unambiguous) determination of the poly- 
nomial Pn-i(xs) we get the following 

Theorem 1. The weight function of the finite- dimensional set of Askey-Wilson polynomials for 
q a primitive root of unity (1.12) and the constraints upon the parameters (2.8), (2.9) has the form 

hN-i{rq' -r-^q-') 



where s — 0,1, . . . , N — 1, and 



^ {9q-^arq^,ar-'q--q)n ^n_ (212) 
^ {ab, ac, ad; q)^ 



This form of the weight function is not convenient for the analysis of its positivity since a truncated 
4933 series (2.12) depending on the discrete argument s is involved into the definition. In the next 
section we shall derive a more explicit form of Ws- 



3. Difference equation and absence of self-duality 

In this section we derive a difference equation for generic AWP for q a root of unity and show 
that a complete self-duality property (a permutational symmetry between the parametrization of 
polynomial's argument and its degree) is absent. 

In order to derive this difference equation we use contiguous relations for the terminating basic 
hypergeometric series i^pz- Recall that if one denotes 

where one of the parameters ^, 77, 7, S is equal to q~'^ , M a positive integer, and efh — ^r]^5q, then 
one has the following relation [2, 9]: 

Bi4>{^q~\ m) + 530(e, V) + B2<t>{^q, r^q'^) = 0, (3.1) 

where 



B, = n{l-r]){iq-r^){i-e){i-f){i-h), (3.2) 

B2 = ^{l-m-m){v-e){v-f){v-h), (3.3) 

^3 = iri{i-r^q){iq-ri){i-v){l-l){l-5)-B,-B2. (3.4) 

In our case one should choose the parametrization 

^ = arq^ , rj — ar~^q~^ , 7 = 5 — gq^~^, e — ah, f — ac, h — ad, (3-5) 
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which gives Pn{xs) = Dn(f){arq^,ar Substituting (3.5) into (3.1) we arrive at the difference 

equation for AWP with the argument Xs given by (2.6): 

A,PniXs-l) + CsPniXs+l) - (As + Cs)Pn{Xs) = A„P„(2;,), (3.6) 

where 

^ [q- - - gq--^) (3.7) 

and 

gq~^{l — ra~^q^){l — rb^^q^){l — rc~^q^){l — rd~^q^) 

' ^ (1 -r2g2s-l)(l _^2^2s) 

(1 — arg'')(l — brq''){l — crg'')(l — drq"") 
' ^ (l-r2g2s)(i_^2^2s+i) ■ 

Let us note that after the (complex) shift of the variable s such that ra~^q^ — > q'^ and a change 
of notation for parameters, the equation (3.6) coincides with the three-term recurrence relation for 
AWP. This means that the difference equation (3.6) may be considered as a recurrence relation for 
a special class of associated AWP. The general associated AWP are determined by (1.1), (1.2) and 
the coefficients m„, obtained from (1.6) after the replacement of g" by //g", where // is an arbitrary 
complex parameter [11]. 

But the recurrence coefficients As and Cs do not have zeros for any s = 0, 1, . . . , so that one can 
consider orthogonal polynomials of unlimited degree. The orthogonality measure of these polynomials 
will be continuous, i.e. the parameter on the right hand side of (3.6) will not be quantized. 
Moreover, there exist regions of the parameters a, 6, c, d such that the coefficients Un are positive and 
hn are real. For instance, such a situation takes place if we set 

a = g^ d^q\ b^-q^, c ^ -q\ g = e^-/^ 

where a, (3, 7, S are real nonzero parameters, choose N even, and impose the restrictions 

\a\, I7I, <S 1, a <0, a < P < —a, 7 > —a, S > —a. (3-11) 

Under these conditions the difference equation (3.6) becomes Hermitian. From this fact follows the 
existence of g-orthogonal polynomials (a special subclass of the associated AWP) for q a root of unity 
with the continuous positive measure which are different from the sieved polynomials. Consideration 
of these polynomials lies beyond the scope of present work. 

Before calculating the weight function, let us discuss a self-duality of the AWP [5, 14]. For 
simplicity let us limit consideration to finite sets of orthogonal polynomials with positive measures 
Pn{s), where n is the degree of polynomials and s is the index numbering discrete values of the 
argument Xg, n,s — 0,1, . . . , N — 1 > 0. Due to the Favard theorem, any set of such polynomials 
is uniquely determined by the three-term recurrence relation. As usual, one has some discrete 
orthogonality relations (1.4) and dual conditions: 

1 = 

n=0 "-n 

(see, e.g., [4]). In the dual picture, P„(s) are considered as eigenvectors with eigenvalues A„ of some 
Hermitian matrix acting upon the (dual) index s. The latter matrix is not necessarily of the Jacobi 
form; in general it is difficult to find explicitly. Moreover, it is not unique (there is a dependence 
on functional parametrization of Xg on s). When this matrix is tridiagonal, one may refer to the 
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(3.8) 
(3.9) 



(3.10) 



(3.12) 



corresponding polynomials as partially self-dual, because then both indices n and s can be considered 
as degrees of orthogonal polynomials. Polynomials Pn(s) arc called completely self-dual if there exists 
a symmetry such that Pn{s) = Ps{n) (this relation requires a special normalization of polynomials, 
see [5, 14]). Here the involution sign means a permutation of parameters of polynomials. In 
the simplest case there is no transformation of parameters at all and the symmetry condition is 
Pn{s) — Ps{n), which means that Pn{s), being considered as a x matrix, equals its transpose. 
For completely self-dual systems, the three-term recurrence relation in n and the difference equation 
in s coincide up to a change of parameters. In the P„(s) = Ps(n) case, these equations just totally 
coincide. For examples of such a situation we refer to the discrete series representations of the 
quantum algebra S0q{3) [17] describing g-ultraspherical polynomials at g"^ = 1. 

In [5, 14] completely self-dual orthogonal polynomials were considered and it was proven that the 
only polynomials possessing this property are the g-Racah polynomials. The recurrence coefficients 
^n, Vn in (1-6) for these polynomials [2] satisfy the conditions .^o = ^a^+i = 0. For the case = 1 
with generic values of parameters only the first factor ^„ creates zeros = ■^Af = 0. I.e., our Pn{xs) do 
not belong to the g-Racah polynomials and hence they are not completely self-dual. The latter can 
also be seen directly from the explicit expression (1.9): for t = rq^ the permutation of n and s cannot 
be compensated by an appropriate normalization of polynomials and permutation of parameters. In 
fact, the values of variable s can be extended from a finite set s = 0,1, . . . , N — 1 to all integers 
without any problem. As a result, one gets the periodicity property Pn{xs+N) = Pnixs) and the 
same periodicity is satisfied by the coefficients Ag,Cs in the difference equation (3.6). Evidently, 
these As, Cg never vanish. AgCs-i 7^ 0, so that there is not automatic reduction of the difference 
equation to the finite-dimensional matrix form in contrast to the g-Racah polynomials case. 

As seen from (1.9), the polynomials P„(xs) become completely self-dual if the condition r = aq^ , 
j integer, is fulfilled. But this equality holds only if one of the products ab, ac, ad is equal to g^, k 
integer, as it follows from (2.5). However, such constraints correspond to the g-Racah polynomials 
and they are forbidden in our case (2.9). To summarize the above considerations: duals of the AWP 
for = 1 are given by the special types of the associated AWP for q^ — 1, i.e., there is no complete 
self-duality. 

Let us return now to the problem of calculation of the weight function Wg- Existence of the 
difference equation (3.6) allows us to find a simpler form of Wg than (2.11). Indeed, multiplying 
equation (3.6) by a function a{s) satisfying the only requirement of periodicity a{s + N) = o"(s), and 
combining it with a similar equation for the polynomials Pm{xg), one can get a bilinear relation: 

Aga{s)[Pn{Xg^l)Pm{Xg) - Pn{Xg)Pm{Xg^i)] 
+Cga{s)[Pn{Xg+i)Pm{Xg) - Pn{Xg)Pm{Xg+l)] 

= (An - Xm)cr{s)Pn{Xs)Pm{Xs)- (3.13) 

Choose a{s) in such a way that 

Ag+ia{s + l)^Cga{s). (3.14) 
Then the equation (3.13) can be rewritten in the form 

Rnm{s + 1) - Rnm{s) = (A„ - X^) CT (s) Pn{x g) Pm{x g) , (3.15) 

where 

Rnm{s) = Aga{s)[Pn{Xg)Pm{Xg--i) - Pn{Xg^i)Pm{Xg)]. (3.16) 

Summing (3.15) from s = Otos = A'" — 1 and using the obvious periodicity property Rnm{s + N) — 
Rnm{s), one gets the orthogonality relation: 

N-l 

(An - A„,) J2 Cr{s)Pn{Xs)Pm{Xs) =0, U ^ TU. (3.17) 
s=0 
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The weight function a{s) is determined from the relation (3.14) uniquely, up to a multiplicative 
factor. Hence one should have Wg = k,(t{s) where k is some constant and w{s) is the weight function 
calculated in (2.11). From (3.14) we find: 

w = w ( / Y {'^-r^l'^'){(^r,br,cr,dr;q)s 
* (1 — r'^)[qr/a, qr/b, qr/c, qr/d; g)^' 

where wq is a normalization constant. Note that the consistency condition a{N) — a{0) imposed 
upon g-difference equation (3.14) for q^ = 1 (see [6] for a general discussion of such conditions) is 
satisfied in a quite nontrivial manner due to (2.7). 

We have found thus the explicit dependence of the weight function Ws on the discrete argument 
s. The remaining value Wq can be evaluated from (2.11) 

^ N-l (1 - ibcr){l - {cdr){l - (M)^)(l - 9/q){l - 9/q'){r - r"^) 
° F{0){l-g^)^{l-bc/q){l-cd/q){l-bd/q){r^ -r-^) ' ^ ' 

where F{0) is fixed by (2.12) at s = 0. The normalization condition 

N-l 

(3.20) 

s=0 

where Ws is given by (3.18), leads to some non-trivial trigonometric identities which we summarize 
into the following theorem. 

Theorem 2. Under the constraints (2.8), (2.9) upon the complex parameters a,b,c,d and for 
q — e^TT^Af/JV^ _ following trigonometric identity takes place 

v^Y^y i'^ - r^(l'^')i(^rM, cr, dr; q)s ,^21) 
s=o U/ {qr/a,qr/b,qr/c,qr/d;q)s 

{ab, ac, ad, be, cd, bd; g)jv-i(l - fl'/?)(l - q/q'^) ' 

where 

, n _ {ab, ac, ad; q)N-i ^ (gq^^, ar, ar^^; g)„ „ 
na, b, c, d) = i; . (3.22) 



Comments. Probably the truncated ^ip^ series (3.22) has a simpler representation through some 
finite products due to the specific properties of the roots of unity and algebraic relation between r and 
a, b, c, d, but we were not able to find it. Since the AWP are symmetric with respect to permutation 
of parameters, the same symmetry should be valid for /(a, b, c, d), but it has not been proved yet to 
be explicit. Since the left hand side of (3.21) is finite for a ^ 0, the function (3.22) should be finite in 
this limit as well. The above identity looks as a root-of-unity analog of a special case of the Watson 
transformation formula for 8</?7 series [9], but the complete coincidence is not yet established. 

The relation (3.21) exists only for q a root of unity, i.e., there are no its analogs for other 
values of q contrary to the g-Racah polynomials' weight function normalization condition holding for 
|g| < l,\q\ > 1, or |g| = 1. As shown below (see also the Appendix), at the bottom of trigonometric 
identities obtained for q'^ = 1 like (3.21), one finds Gauss sums or similar summation formulas. 
The latter play an important role in the number theory. Therefore, the arithmetic aspects of (3.21) 
require a detailed investigation. Note in passing that the most general identity of such type should 
be tied to the biorthogonal rational functions associated with the special io</?9 series [10]. 
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4. Two parameter symmetric Askey- Wilson polynomials 



In this section we consider a special subclass of the symmetric AWP, and determine some regions of 
parameters for which the measure is positive. 

Consider AWP with a = — c = g", d = —h = . Then in the recurrence relation for AWP only 
nondiagonal terms survive, i.e., the polynomials are symmetric: 

Pn+l{x) + UnPn-l{x) = xPn{x), (4.1) 

where n = 0, 1, . . . , iV — 1, 

4sina;n smujiji + 2a + 2(3 — 2) cosuj{n + 2a — 1) cosu{n + 2/? — 1) , , 

~ sm2uj{n + a + (3 - 3/2) sin 2a;(n + a + /3 - 1/2) ^ ' ' 

and u = ttM/N. 

One can show that if M = 1 and is odd, the recurrence coefficients Un are positive, Un > 0, for 
n — 1,2, . . . , N — 1, under the following constraints upon the parameters a and (3: 

-1/4 <a< 1/4, 1/4 < (3 < 3/4, (4.3) 

or 

1/4 < a < 3/4, 3/4 < /3 < 5/4. (4.4) 

These restrictions, or their partners obtained after permuting the parameters a and (3, provide only 
particular examples of positive Un, the problem of finding all such cases (including the cases of M 7^ 1 
and even N) is complicated and is not considered here. 

The difference equation for the taken symmetric AWP has the form (3.6), where 

^ _ sin2a;(s-a + l/4) sin2cu(s-/3+l/4) 

' ~ sin2u;(s + 1/4) sm2u{s - 1/4) ' ^ ' 

^ _ sin2a;(s + a + 1/4) sin 2a;(s + + 1/4) 

' ~ sin2u;(s + l/4) sin2u;(s + 3/4) ' ^'^ 

Xn — — 4sina;n sina;(n + 2a + 2(3 — 1), 

Xs = 2cos2u;(s + 1/4). 

Under the conditions (4.3) or (4.4), this difference equation can be transformed by a gauge transfor- 
mation Pn{xs) = x{^)4'n{s) with an appropriate periodic function x(s) to the hermitian form 

Us+iMs + 1) + U,iPn{s - 1) - {A, + Cs)Ms) = KMs), (4.7) 

where Us = y/A^U^ is a real function because AgCg-i > for all integer s. Hence (4.7) provides 
a nontrivial example of the Hermitian Jacobi matrix with periodic coefficients whose cigenf unctions 
and eigenvalues are known explicitly. Note that in this case, the difference equation for symmetric 
polynomials is determined by the tridiagonal Jacobi matrix (cf. [17]). 
The weight function for the corresponding AWP has the form 

„. ^ , . /.X A sin2w(/c + a-3/4)sin2tu(A; + /?-3/4) , 

Rsm2u;(s+1A)Y\ }- '-^ }- ^ Vr, 4.8 

^ ' ' ^}^sm2uj{k-a + l/A)sm2uj{k- P + l/A)' ^ ' 

where s = 0,l,...,A^ — 1, and i? is a normalization constant. 

One can check that > for the intervals (4.3) and (4.4) as it should be (by the Favard theorem) 
for the case of positive recurrence coefficients As in the generic case, the considered symmetric 
polynomials do not obey complete self-duality. 
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5. Finite orthogonality relation for continuous g-Jacobi poly- 
nomials 

In this section wc study a finite set of the continuous g-Jacobi polynomials [3, 9, 12] for g a root of 
unity. These polynomials appear from the AWP with arbitrary parameters a and c and the choice 

b = aq^/^ d = cgV2 (5.1) 

(we use a slightly different parametrization than in [3, 9, 12]). Moreover, we will assume that M in 
(1.12) is odd. Then from (2.5) we see that En — 0, hence the zeros of the A'"th continuous g-Jacobi 
polynomials are 

Xs = 2 cos7r^(s + 1/2) = g(^+V2)/2 + ^-(.+i/2)/2^ s = 0, 1, . . . , TV - 1. (5.2) 

We could set M = 1 in (5.2). However this parametrization is more convenient for further cal- 
culations. The zeros (5.2) are simple, real, and they coincide with zeros of the monic Chebyshev 
polynomials of the first kind Tjv(x) = 2 cos NO. 

In order to simplify the expression for Pn_i[xs), we use another form of the continuous Jacobi 
polynomials due to Rahman [3, 9] 



where p = q^l'^ = e^''^^^ and 



" p'^/'^[—acp'^]p)n 



The formula (5.3) is vahd for all n = 0, 1, . . . , A?" — 1. However, for n = — 1, a remarkable 
simplification takes place. Since M is odd we have p^~^ — —p, hence (5.3) is reduced to the 
balanced 3(^2 function 

p / X n (p-',p'+\acp-^ \ 

But this sum can be calculated exactly using Jackson's extension of the Pfaff-Saalschutz formula [9] 
and adapting it to roots of unity. After this, we have 

Pn-.{Xs) - Dn., \-j (,^V2;p)^(epl/2;p), " (^'S) 

The Pn{xs) can be easily calculated because Pn{x) — 2cos9N, x — 2cos^. Namely, we have 

^^^^^^ = sin(7rM(. + l/2)/Ar)- ^^"^^ 
From (1.5) we find thus the weight function of the continuous g- Jacobi polynomials 

^ R^M^Mis + l/2)/A^) ^7^Tll''C'^y3^V^^ ^ = 0,l,...,Ar-l, (5.8) 
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where the normahzation constant is 

Rn = ip ' — 7 T 7 7 • (5.9) 

(-acp;p)iv-i(-p;p)iv-i 

In general, the parameters a, c are complex and the weight function (5.8) is also complex. Nev- 
ertheless, there are important cases when the weight function becomes positive. Indeed, put M — 1 
and choose the following parametrization [3] 

a = p"+V2, (5.10) 

which is quite natural for it leads to the ordinary Jacobi polynomials P^'^\x) when N oo. Using 
this parametrization and passing to the trigonometric notations, we can state the following theorem. 

Theorem 3. The weight function of the finite- dimensional q-analogs of the Jacobi polynomials 
(5.3) for q — exp{2Tri/N) has the form 

„ , , , , sina;(A; + a) cosa;(/c + /3) 

where ui = n/2N, s = 0, 1, . . . , — 1, and the normalization constant 

Rj^ = W siJ^^(^ + /^)cos^(fe + Q) ^2) 
k=i cosu}ksinu}{k + a + /? + 1) 

For 

-Ka, P<1, (5.13) 
the weight function (5.11) is positive, Ws > 0. 

Proof of the positivity of Wg is straightforward. In this case, the recurrence coefficients Un are 
positive too and the coefficients 6„ are real. So the condition (5.13) leads to the class of ordinary 
orthogonal polynomials [8] having the positive weight function on the real spectral set. The particular 
choice a = P corresponds to the g-ultraspherical polynomials considered (for the case of g a root of 
unity) in [17]. There are also nine special cases when the weight function (5.11) is reduced to simpler 
expressions. These cases arise when a,P — —1/2,0, 1/2. For example, when a — P — —1/2 we get 
finite set of Chebyshev polynomials of the ffist kind with the constant weight function Wg = 1/N; for 
a = f3 = we obtain Legendre g-polynomials with = sincj sinu;(2s + 1); for a = —1/2,(3 = 1/2 
we have the weight function Ws = iV""*^ cota;/2cos^ci;(s + 1/2), etc. 

Note that the normalization condition Ws — 1 can be considered as a nontrivial summation analog 
of the standard beta- integral, because in the limit g — > 1 (i.e., — > oo) we get the ordinary Jacobi 
polynomials with continuous weight function w{x) = (1 — x)°'{l + x)'^. Normalization of this function 
leads to the ordinary beta-integral, g- Extension of the beta integral was found in [3] when q is inside 
the unit disk, \q\ < 1. 

Consider some limiting cases. When c ^ 0, we get the continuous g-Laguerre polynomials [12] 
with the weight function 

^ + l/2)/N)a--p--y^ /iy- . (6.14) 

{-P;P)n-i [a-^p^l^-p)s 



When a = p^^^ , we get an interesting trigonometric identity 

N-l 

sin7r^(s + 1/2) p-(-+i)/2 ^ _^pi/2_ 
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If in (5.14) we take the limit a — > then we get the weight function for g-Hermite polynomials 

/ , sm(TrM(s + l/2)/N). (5.15) 

(-p;p)jv-i 

6. Finite orthogonality for big g-Jacobi polynomials 

Big g-Jacobi polynomials can be considered as some limiting subcase of the AWP [9]. However, it 
is more convenient to work directly with these polynomials themselves. Their recurrence coefficients 
are 

where 

^ +i (l-g")(l-&g")(l-a&c-V) 

^ (l-a6g2n)(l_ai,^2n+l) ' 

(1 - ag"+^)(l - cg"+^)(l - abg"+^) 

(l-a6g2n+i)(i_a6g2n+2) ■ ) 

The explicit form of the big g-Jacobi polynomials is 

V cq J 

where 

(aq; q)n{cq; q)n 



(abq^+^;q)n 

For g fixed by (1.12), the A'"th polynomial has the form 

In what follows, we set c = 1. Then the zeros of the polynomial Pn{x) coincide with the roots of 
unity and lie at the vertices of the regular polygon 

x, = g^ s = l,2,...,7V. (6.4) 

For Pn^i{xs) we have 

Pn-i{xs) = Dn_i 2<fi ;q,q] ■ (6.5) 

\ aq J 

This expression can be evaluated using the root-of-unity analog of the Chu-Vandermonde formula 
[9] yielding 

Pn-i{xs) = D^^,{abf-^^^^^^^, (6.6) 

{aq;q)N-s 

where 

(l-a^)(l-a&g-^) 
^^-^ = '^(l-a)(l-(a6)^) ■ ^'-^^ 



Hence for the weight function, using (1.5), we have 

(l-a^)(l-a6g)(6;g),g^ 



aq(b-l)(l-a^b^)(a-^;q) 
11 



-, s = l,2,...,iV. (6.^ 



It can be verified that this weight function is normalized, J2^=i Ws = 1- 

There arc two interesting limiting cases: when a — > oo we get the g-Meixner polynomials [12] 
with the weight function 

Ws = b'-''il-b)-\b;q),q^. 
When 6 — > we get the big g-Laguerre polynomials [12] with the weight function 

Ws = {qa)~^{a^ — V 



It should be noted however that in both these cases we obtain only special (one-parameter) type of 
the corresponding polynomials due to the condition c—1. 



7. Finite orthogonality for alternative g-Jacobi polynomials 

So far, wc were considering AWP with the base q. In this section we consider a special class of the 
AWP with the base p = q^l'^ = e^'^^^l^ ^ where M is odd. These polynomials are orthogonal on some 
points of the real interval and can be interpreted as alternative g-Jacobi polynomials. 

Let us set c = —d = 1 and keep a and b as arbitrary complex parameters for AWP. Then the nth 
Askey- Wilson polynomial (with the base p) is written as (cf. (1.9)) 

p/x n (p~'',-abp''-\az,a/z \ 
\ ao, a, —a J 

where x = z -\- and in the expression for normalization constant Dn (1-10) wc should substitute 
p instead of q. These polynomials do not belong to the family of finite polynomials considered in 
Section 2, but rather they are subcases of the g-Racah polynomials considered in [2]. Still, as far as 
the authors understand, the specific property of these polynomials to have positive weight functions 
for some particular choices of M and has not been discussed in the literature (see also [16, 17]). 
Let n = A^ + 1. Then taking into account that = — 1 we get 

TD r \ n (p-^-\az,a/z \ 

\ a,-a J 

where 

i^iv+i = a-"^-^(l-a^)(l-a^^). 
In order for polynomials P„(a;), < n < A" + 1, to be well defined we impose the constraints 

ab,a,b^p^, A; = 0, 1, . . . , 2A^ - 1. (7.3) 

In particular, in this case Un does not vanish between two zeros at n = and n — N -\- 1. 
The expression (7.2) can be simplified using the Jackson summation formula [9] 

Pm+i{x) = {z;p)n+i{z'';p)n+i = z-^'-'iz' - l){z''' - 1). (7.4) 

Prom (7.4) we find zeros Xg of the polynomial P^v+i (x) 

Xs=p' +p~' = 2cosnMs/N, s = 0,1,2, N. (7.5) 

It is seen from (7.5) that all the zeros Xg are real and simple. 
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Prom (7.4) one can derive 

P'^^,{xs)^2N{-ir/Xs. (7.6) 
where Xs is a "characteristic function" defined as follows 



Xs 



ri, s = i,2,...,Ar-i 

11/2, s = 0,7V. 



In order to calculate Pn{xs) note that AWP are invariant under any permutation of the parameters 
a, 6, c, d. Hence we can write 

' ab '^'^J ' 

where 

^ ^ (-l;p)iv(Q;p)jv(&;p)iv 
{abp-^;p)N 

Again the expression (7.7) can be simplified using the Jackson summation formula: 

Pn{xs) = Dn—, r-fT—\ P) ' s = 0, 1, 2, . . . , AT. 7.8) 

Combining all these results together we can state the following theorem. 

Theorem 4. The weight function of the finite set of alternative q-Jacobi polynomials (7.1) with 
p — exp{TriM/N), {N, M) — 1, M odd, and restrictions upon parameters (7.3), has the form 

hN . {a;p)s{b;p)s{~p/aby 

Pn[xs)Pn+i{Xs) {p/a;P)s{p/b;p)s 

where s = 0,1, . . . , N and 

^ ^ {p;p)N{-a;p)N{-b;p)N 
2N{-ab;p)N 

The normalization condition X^fLo'^s = 1 in this case leads to some interesting trigonometric 
identities. Indeed, taking the limit a — >^ 0, 6 ^ we obtain the weight function 

= Xs (-1) P ■ (7.10) 



The normalization condition leads to the celebrated Gauss sum [7] (see also the Appendix) 

E - {-p;p)n-. = V]vp^(^-)/^ (^) , (7.11) 



where (-^) is the Jacobi symbol. For M — piP2 • ■ -Pk ~ the prime number decomposition of M - 
this symbol is defined as the product 

m) \pi) \pk)' 

where (^) is the Legendre symbol. The latter is equal to +1 if A" is a quadratic residue (i.e., if there 

exists an integer x such that x'^ — N (mod pj)), to —1 if A" is not a quadratic residue, and to if pj 
divides A". In the simplest case M — 1, one has (y) = 1. 
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In another limiting case, setting a — p^/^ and taking the hmit 6 — > 0, we get the weight function 

Ws = A^Xsp-'"^\ (7.12) 

where 

. {P]P)n{-P^^^]P)n {-P^'^]P)n /„ion 
^iV (-p;p).-.- ^'-^'^ 

Using the formulas (8.16) and (8.19) from the Appendix, one can simplify the expression (7.13) for 

--/l^"'1f)- (-^) 

Hence we have the following trigonometric Gauss-like sum 

IL^P-''" - \{l - (-0"") + /f exp(-i.M/4) (^) . (7.15) 

For M — 1 this identity was derived in [13]. The M ^ 1 cases were probably discussed in the 
literature as well. 

At the final stages of preparation of this paper, the authors have discovered that in 1968 M. 
W. Wilson had considered a discrete analog of the Legendre polynomials [18], which is a particular 
subcase of the above alternative g-Jacobi polynomials. Indeed, let us set in (7.9) a — —b — p^"*"^ 
and take the limit e — > 0. As a result, one finds the weight ftmction wq = wn = and Wg — 
tan(7rM/2A^) sin(7rMs/A^) for s = 1, 2, . . . , — 1. For M = 1, it coincides with the weight function 
of [18] up to a normalization constant. Although some important characteristics of the corresponding 
polynomials (e.g., the nice approximation properties) have been found in [18] the description was not 
complete. In particular, the explicit form of the recurrence coefficients, which has not been derived, 
is 

_ sin2 7rn/2A^cos7r(n - l)/2A^cos7r(n + l)/2Ar , _ 

~ sin7r(2n-l)/2A^sin7r(2n+l)/2Ar ' ~ ^ 

Only the first three coefficients ^1,2,3 were found in [18] and they coincide with (7.16). Note that the 
corresponding polynomials are g'-analogs of the Legendre polynomials different from those mentioned 
above which were considered in [17]. The latter are simpler but it is natural to expect that they 
provide a similar approximation of the Legendre polynomials for N ^ 00. 

Returning to the general weight ftmction (7.9), we note that when a = p", b = —p^, and < 
a, /9 < 1, the recurrence coefficients m„ become positive leading to the positivity of the corresponding 
weight function w^, s = 0, 1, . . . , A^. Verification of this statement is elementary. In the limit p — > 1 
we recover the ordinary Jacobi polynomials P^'^'^''^\x). We have thus a subclass of the "classical" 
orthogonal polynomials (in the sense of positivity of the weight function on a real interval) among 
the alternative g- Jacobi polynomials in the same way as for the case of the continuous g'- Jacobi 
polynomials. 

Although we did not give a complete classification of the regions of complex parameters a,b,c,d 
and q in the recurrence coefficients (1.6) leading to positive measures for the corresponding (associ- 
ated) Askey- Wilson polynomials, the present work provides some necessary tools for the resolution 
of this important problem. 

The authors are deeply indebted to R. Askey and G. Gasper for helpful remarks to the preliminary 
version of this paper. 
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8. Appendix 

In this appendix we present some useful formulas concerning g-calculus when g is a root of unity. 
Some of them are well known, nevertheless they are included here because the systematic treatment 
of such formulas is absent in the literature. We consider here only standard g-hypergeometric series, 
although there are similar relations for the bilateral series as well. 

Assume that g is a primitive root of unity, i.e., q = e^'^*^^/^, where {N,M) = 1. Denote also 
p = q^l"^ = e'^*^/^. The first set of formulas is connected to the g-shifted factorial defined in the 
standard manner 

(a; g)„ = (1 - a)(l - ag) . . . (1 - aq^-'), (a; q)o = 1. (8.1) 

Obviously {a;q)n+N — ?)n(Q; so it is sufficient to study only the case < n < N. The 
following identities are easy to derive 

(a;g)^ = l-a^ (8.2) 
(in particular (a; g) jv = if a is a primitive A^th root of unity) and 

1 - 

{a]q)N-i^- TT' iia^q, (8.3) 

1 — aq 

{q;q)N-i = N. (8.4) 

In general, there are no simple formulas for {a;p)N and {a;p)N_i like (8.2) and (8.3). In our case the 
following relations for inversion of the parameter are valid. Assume that M is odd. Then 

(a;p)^ = aV^"^+^)/'(-a-V;p)iv, (8.5) 
{a;p)N-i = a^-y-^(^+^)/^(-a-y;p)^_i. (8.6) 

If M is even (hence A^ is odd), then 

{a;p)N = -a^{a~^p;p)N, (8.7) 
{a;p)N-i = pa^"^(a"V;p)Af-i- (8-8) 

The next set of formulas is connected to finite sums involving the g-shifted factorial. The simplest 
one is a special case of the g-binomial theorem which can be written as 

N-s ( s. „\ 1 _ 

where s = 1, 2, . . . , A" — 1. In (8.9) and below we use the standard notation for g-hypergeometric 
functions although we always assume that only the first A" terms are present in the corresponding 
series. 

The Chu-Vandermonde formula [9] for q a root of unity looks like this: 

where a,c are arbirary complex parameters, ^ 1. In particular, for s = 1 we get a pretty 
summation formula 

^Ma;g), (a^-c^)(l-cg-^) 

h (c;?)^ (l-c^)(a-cg-i) • ^ - > 
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Similarly, Jackson's extension of the Pfaff-Saalschutz sum [9] in this special case is written as 

{a;q)k{b;q),q' _ (g^ - c^)(6^ - c^)(g - c)(a6g - c) 
h (c; q)k{abc-^q'; q)k (1 - c^){a^h^ - c^){qa - c){qh - c)' > 

The Dixon ^-summation formula [9] takes the form 

{-qa-q),{h-q), , {h - a^){l + a^){h^ - a^) 

A special case of the Singh quadratic transformation yields, for — 1, the summation formula 

\a^,-ah ^q,-ab V / {-a;p)N-i{b/a;p)N-i 

where we assume that M is odd so that = —1. If one takes the limit a — >■ 0, 6 — >■ 0, 6/a — >■ then 
we get an important summation formula 

Af-l 

E(-i)V/' = (-p;p)iv-i, (8.15) 

k=Q 

going back to the works of Gauss in the beginning of nineteenth century. The Gauss sum (8.15) 
plays a prominent role in some problems of the number theory [7]. The right hand side of (8.15) can 
be rewritten as follows [7] 

(-p;p)^-l=(^)v^e-^(^-^/^ (8.16) 
where (^) is the Jacobi symbol. Analogously, considering the Gauss sum 

2N-1 

Y.i-irp-''^' = i-p'^';p'%N-i, (8.17) 

it=0 

one can get the identity 

i-p'^'-y%N-i = (^) V^e-^(^^-^)/l (8.18) 
Prom this relation and the obvious equality 

{-p;p)n-i {-p"';p)n = {-p"'-y%N-i 

follows another interesting identity 

where the following properties of the Jacobi symbol were used 



16 



References 

[1] W. Al-Salam, W. R. AUaway, and R. Askcy, Sieved ultraspherical polynomials, Trans. Amer. 
Math. Soc. 284 (1984), 39-55. 

[2] R. Askey and J. Wilson, A set of orthogonal polynomials that generalize the Racah coefficients 
or6-j symbols, SIAM J. Math. Anal 10 (1979), 1008-1016. 

[3] R. Askey and J. Wilson, Some basic hypergeometric orthogonal polynomials that generalize Jacobi 
polynomials, Mem. Amer. Math. Soc. 54 (1985), 1-55. 

[4] F. V. Atkinson, "Discrete and continuous boundary problems" in Mathematics in Science and 
Engineering, 8, Academic Press, New York, 1964. 

[5] E. Bannai and T. Ito, Algebraic Combinatorics I. Association Schemes, Benjamin/Cummings, 
Menlo Park, 1984. 

[6] R. D. Carmichael, The general theory of linear q-difference equations, Amer. J. Math. 34 (1912), 
147-168. 

[7] K. Chandrasekharan, Introduction to Analytic Number Theory, Grundlehren Math. Wiss. 148, 
Springer- Verlag, New York, 1968. 

[8] T. S. Chihara, An Introduction to Orthogonal Polynomials, Math. Appl. 13, Gordon and Breach, 
New York, 1978. 

[9] G. Gasper and M. Rahman, Basic Hypergeometric Series, Encyclopedia Math. Appl. 35, Cam- 
bridge University Press, Cambridge, 1990. 

[10] D. Gupta and D. Masson, Contiguous relations, continued fractions and orthogonality, to appear 
in Trans. Amer. Math. Soc. 

[11] M. E. H. Ismail and M. Rahman, The associated Askey-Wilson polynomials. Trans. Amer. Math. 
Soc. 328 (1991), 201-237. 

[12] R. Koekoek and R. F. Swarttouw, The Askey scheme of hypergeometric orthogonal polynomials 
and its q-analog. Delft University Technical Report 94-05, 1994. 

[13] N. M. Korobov, Exponential Sums and their Applications, Math. Appl. (Soviet Ser.) 80, Kluwer, 
Dordrecht, 1992. 

[14] D. A. Leonard, Orthogonal polynomials, duality, and association schemes, SIAM J. Math. Anal. 
13 (1982), 656-663. 

[15] A. F. Nikiforov, S. K. Suslov and V. B. Uvarov, Classical Orthogonal Polynomials of Discrete 
Variable, Springer Ser. Comput. Phys., Springer- Verlag, Berlin, 1991. 

[16] V. Spiridonov and A. Zhedanov, Discrete Darboux transformations, the discrete time Toda lat- 
tice, and the Askey-Wilson polynomials. Methods Appl. Anal. 2 (1995), 369-398. 

[17] V. Spiridonov and A. Zhedanov, q -Ultraspherical polynomials for q a root of unity, Lett. Math. 
Phys. 37 (1996), 173-180. 



17 



[18] M. W. Wilson, On a new discrete analogue of the Legendre polynomials, SIAM J. Math. Anal. 
3 (1972), 157-169. 

V. Spiridonov: Centre de Recherches Mathematiques, Universite de Montreal, CP. 6128, succ. 
Centre-ville, Montreal (Quebec) H3C 3J7, Canada; current address: Laboratory of Theoretical 
Physics, JINR, Dubna 141980 Moscow region, Russia 

A. Zhedanov: Donetsk Institute for Physics and Technology, Donetsk 340114, Ukraine 



18 



